String order and symmetries in quantum spin lattices 
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We show that the existence of string order in a given quantum state is intimately related to the 
presence of a local symmetry by proving that both concepts are equivalent within the framework of 
finitely correlated states. Once this connection is established, we provide a complete characterization 
of local symmetries in these states. The results allow to understand in a straightforward way many 
of the properties of string order parameters, like their robustness/fragility under perturbations and 
their typical disappearance beyond strictly one-dimensional lattices. We propose and discuss an 
alternative definition, ideally suited for detecting phase transitions, and generalizations to two and 
more spatial dimensions. 
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Order parameters play a crucial role in describing the 
different phases of matter. However, there exist some 
phases which, despite displaying very intriguing features, 
are not amenable of such a description. In some cases it 
is nevertheless possible to introduce more sophisticated 
quantities that are able to characterize those phases. A 
paradigmatic example is given by the string order pa- 
rameter (SOP), which reveals the appearance of a hidden 
order (so-called string order, SO) in certain spin systems 
[EH Hi- This quantity can be expressed as an expecta- 
tion value of some non-local operator, and the appear- 
ance of SO is highlighted by a non-vanishing value in the 
thermodynamical limit. Despite its importance, we still 
do not have a systematic characterization of its proper- 
ties. It is not clear under which conditions SO appears 
in a ID system, for which kind of small perturbations re- 
specting the gap it is robust [J] , or when it can be used to 
detect a quantum phase transition. Apart from that, it 
seems that the SO looses some of its desirable properties 
beyond strictly one-dimensional systems [H, 0] • 

In this work we clarify all those questions for finitely 
correlated states (FCS) [f|, i.e., matrix product states 
(MPS) i on infinite chains. The relevance of these 
states relies on the fact that every quantum state of a 
finite system has an exact MPS representation [H, 0] and 
that ground states of ID short-range interactions can be 
efficiently approximated within this class In this 

framework we will show that the appearance of SO is 
intimately related to the existence of symmetries, which 
explains how it can be used to detect quantum phase 
transitions. We propose another parameter which better 
recognizes this appearance, since it does not have some 
of the shortcomings of the SOP. We also provide a nat- 
ural generalization of SO to higher dimensional lattices 
(membrane order) which retains all the desired proper- 
ties. Finally, we give several examples displaying a large 
variety of phenomena. 

ID Chains: String order. We will mostly consider 



infinite chains of identical spin-S 1 particles in a transla- 
tionally invariant state '5. We say that the state W has 
string order if there exist a local unitary 1 and local 
operators x, y (which can be taken hermitian) such that 
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Jim \S N (x,y,u,V)\ > 0, 



S N (x,y,v,,^f) 



i^\x <£> u 



(1) 

y|*>- (2) 



Later on we will introduce alternative quantities which 
will extend this definition [l8j |. 

Matrix product states of L spins in a translationally 
invariant state have the form 



nx...ni J =—S 



tr(A ni . ..A nL ) \ni 



(3) 



where the A's are D x D matrices. We will write for 
the thermodynamic limit meaning that L — > oo is taken 
after the expectation value. In this limit the states are 
known as finitely correlated states (FCS)0 — the subject 
of our studies. Most of the properties of these states are 
encoded in a linear map defined as 



S(X) = £ A n XAl 



The matrices A can always be chosen such that 
£(!) = !, £*(A)=A, 



(4) 



(5) 



where £* denotes the map which is obtained by inter- 
changing A n <-> Al in g]), A > and tr(A) = 1. Thus, 
£ {£*) is a unital (trace-preserving) completely positive 
map, i.e., a quantum channel, and has an eigenvalue equal 
to 1 . A FCS is a pure state iff A > and £ has only one 
eigenvalue of modulus one. Since there is a unique de- 
composition [6| of mixed FCS into pure ones and a mixed 
FCS has SO iff one of its pure components has it, we will 
in the following restrict to pure FCS. Note that all of 



2 



them are unique ground states of gapped finite-range in- 
teractions. 

String order and finitely correlated states. For 

any unitary u the SOP © of a FCS is most easily ex- 
pressed by introducing a map 

£ U (X) := Y,(n'Mn)A n XAl =^e^AjXA], (6) 



where A, = £ n <j|n)A„, and u = V ■ e ie i Then | 



S N (x,y > u,9 00 )=ix[A£ a .e^Sy{l)], 



(7) 



where E x , y are defined analogous to £ u in Eq.©. 

The following Lemma studies the spectral radius p of 
£ u which is crucial for Eq. ([7j) due to the limit N — > oo. 

Lemma 1 p(£ u ) < 1 with equality iff there exists a uni- 
tary V and 9 G [0, 2ir) such that 



V^Aj = e^-^AjVl (8) 
£ u has at most one eigenvalue of modulus 1. 



Proof: Let us consider an eigenvector, V, of £ u with 
eigenvalue A, i.e. £ U (V) — XV. Multiplying from the 
right by AV^ and taking traces, we obtain 



|A|tr(FAFt) = 



^2e ie nr(AjVA]AV^ 



< 



^tr^ItAl^t) 
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tr^AW), 



1/2 



J^tr^tvAVUj-) 



1/2 



(9) 



where we have used Cauchy-Schwarz inequality and ([5]). 
Since A > 0, tr(VAW) > and thus |A| < 1 as stated. 
Now, if condition ((5} is fulfilled one can readily see that 
e 10 is an eigenvalue of £ u by using Eq. ©, and thus 
p(£ u ) = 1- On the contrary, if |A| = f, then the 
inequality in ([9]) has to become an equality, and thus 
ae i9 i k 1 / 2 V' [ Aj = A 1 / 2 AjV^. Multiplying by the adjoint 
expression, taking traces, summing in j, and using again 
([5]) one obtains that \a\ — 1, i.e. a — e~ l9 . Since A is 
invertible we obtain ijHJ). This also implies that 



8(V*V) = V^^A^A^V = V^V, 



(10) 



where we have used ([5]). Since 1 is the only fixed point of 
£, we get V'V = 1. Moreover, suppose that £ u has two 
eigenvectors, V, V with eigenvalues e l6) , e , respectively. 
Then, using ([9]) we have 



£(V^V) = ^I^vVlt = eW-VviV, (11) 



such that the same argument gives V = V and 9 = 9' 
□ . 

Now we can specify the conditions required for SO. 
First, p{£ u ) = 1 since otherwise Sn will decay exponen- 
tially with N. Using Lemma 1 we know that the eigen- 
value A of magnitude 1 is not degenerate, so let us denote 
by V and Y the corresponding right and left eigenvectors, 
i.e. £ U (V) = XV, £*(Y) = XY, where £* is again given 
by expression ((6]) but interchanging A n A^. We have 
\im N ^ 00 S N {x,y,u,^ 00 ) = tx\Y £ v {l)]tx[A£j{V)}. Writ- 
ing £*(Y)V = XYV and using §8§ we arrive to the con- 
clusion that Y = AV^ . Thus, the conditions for the 
SOP not to vanish arc: (i) p(£ u ) = 1; (ii) tx[AV^£ y (l)}, 
tx[A£ x (V)] ^ 0. 

We may ask ourselves if the condition p(£ u ) — 1 is 
sufficient to have SO, i.e. if there are always two opera- 
tors x and y such that the other conditions are fulfilled. 
To answer this question we notice that tr[AV 1 '£ y (]L)]* = 



tr[A£,(V)], where 



uy and u = J2 1 



'mi- 



Thus, we can always choose y = x'u so that condition (ii) 
above is simplified to tr[A£' a! (V')] ^ 0. It is clear that it 
suffices that tr(V AAnA^) ^ for some n, m since then 
we can simply choose x = |n)(m|. We conclude that 

Theorem 1 For a pure FCS there exists SO iff there 
exist a unitary u =/= t, V , and n, m such that 



£u(V) = V, tr(VAA n Al) £ 0. 



(12) 



Now, we show that the second condition can be 
dropped in two situations. First, if x and y in ([2]) are 
products of observables acting on D 2 spins. The rea- 
son is that the set Sd ■— span{A ni . . . A nn A\ rlD . . . A\ ni } 
spans the set of D x D matrices, so that there is al- 
ways a linear combination X of these matrices for which 
tx(VAX) 0. To see that this set is complete note first 
that S* m _i C S m since Y, nr> A ni ■ ■ ■ A nD A\ D . . . A^ mi 



A , . . . A nD _ 1 AJ nD ... Al ni . This inclusion must be 
strict unless m > D 2 since S v 



5m- 1 = S m implies S m = 
S m +i and for a sufficient large N the set {A ni . . . A nN } 
must span the entire space of matrices [f| . Another situ- 
ation is the one in which there exists a continuous group 
of unitaries V fulfilling the first condition, i.e., we can 
parametrize V = e i ^ >H . Then, we can always choose 4> 
sufficiently small such that tr(VA) ^ and this suffices 
since we have ([3]). 

Symmetries in finitely correlated states. We say 

that a state has a local symmetry if there is a unitary 
u ^ 1 such that 



u <g> 



I*) 



This formally means that for every iV-site reduced den- 
sity operator g we have u® N gu^ ® N = g. 

For FCS the condition p(£ u ) — 1 is not only equivalent 
to having SO — as we saw in the previous section — it is 
also equivalent to the presence of a local symmetry: 
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Theorem 2 A pure FCS has a local symmetry u iff 

p{£u) = 1- 

Proof: If p(£ u ) = 1 the result is a direct consequence of 
Lemma 1 which implies that V^AV = A. The converse 
follows from the fact that 



1 

d 2 



< tr(£ 2 ) 



tr[gu® N QU^ N ] 



tr 



L(£ U ®£ U ,) N {R) 



for some L, R which are independent of N. □ 

This theorem, together with Lemma 1 provides a com- 
plete characterization of FCS with local symmetries. We 
note that an analogous statement can be found (though 
without proof) in Ref.[g|. 

Equivalent criteria for the existence of a local symme- 
try can be given in terms of the isometry B := J^j |i)A? 
as well as for the D 2 x D 2 matrix E := J2j A? ® Aj • ^ n au 
cases V and u are elements of two unitary representations 
of a symmetry group (different from the identity) : 

• Condition CI: (u ® 1)B = (1 ® V)BV^ . 

• Condition C2: £ is covariant, i.e., for all X, 
£{VXV^) = V£(X)V^. 

• Condition C3: [E, (V <g) V")] = 0. 

CI can be obtained from Lemma 1 by using the spectral 
decomposition of u. C2 and C3 are the same if we use 
(k,l\E\i,j) — (k\£(\i)(j\)\l). It is also clear that Lemma 
1 implies C2. Finally, if C3 is fulfilled, then we have that 
are also Kraus operators of the map £ . Since all 
Kraus decompositions are related by a unitary matrix, 
say u, we have that VAjV^ = ^2, n {n\u\j) A n implying 
the condition of Lemma 1 . 

Now, we can use C3 to derive a criterion for the ex- 
istence of a continuous symmetry where V — e*^ . Ex- 
panding in first order in <f> we obtain 



M(H) := [E, H <g> 1 - 1 <g> H] = 0. 



(13) 



Since M is a linear map on the space of hermitian ma- 
trices, we have that: a pure FCS has a local continuous 
symmetry iff M has a non-trivial kernel. 

Similar criteria can be given for discrete symmetries, 
e.g., Z2 symmetry where the A's are either block diagonal 
or block off-diagonal 

Alternative definitions. One can understand the 
importance of the SOP to detect quantum phase tran- 
sition in terms of its relation to local symmetries. If 
we have a Hamiltonian with certain local symmetry and 
its ground state is unique, then there may be SO. If we 
change the parameters of the Hamiltonian but keeping 
the symmetry until the gap closes, then the ground state 
will be degenerate and the symmetry may be broken [l^ | . 
Thus, the SO may disappear at that point, indicating the 
presence of the transition. Note that due to the possible 
choices of the operators x and y in the definition it 



may happen that the SOP for a particular choice van- 
ishes even if there still is a symmetry. In order to avoid 
this, one may look at the quantity 



Rl(u) 



(*l|u® L |*l 



(14) 



where \^ l) is the ground state of the Hamiltonian acting 
on L sites with periodic boundary conditions (Eq. ([3]) in 
the case of MPS), which is indeed directly related to the 
existence of a symmetry. In fact, i?oo(w) = lim^, Rl(u) 
can only vanish if the gap is closed, and thus it is ideally 
suited to study the presence of transitions. Note also that 
it can be straightforwardly determined from numerical 
algorithms based on MPS [lOfl , 

Example 1: AKLT state. It is instructive to revise 
the appearance of SO in the ground state of the AKLT 
model [lH. For that state we have S = 1, A° = a z /y/3 
and A 1 '^ 1 = y^2/3a± where the sigmas denote Pauli 
matrices. One finds A = 1, and taking u = e z we 
obtain V — o~ z and we can take x = y = S z , so that 
Sifo^Uj^oo) = -4/9 [1]. Note that the AKLT state 
has SU{2) symmetry so that obviously Rl{u) = 1. 

Example 2: Cluster state. We have Q 



.4° 



J_/l 1 

V2 \0 



.4 1 



J_ (0 

V2 V 1 -1 



(15) 



This has the symmetry induced by u — —o~ x . We have 
V = o-y, and one can readily see that tr(V AA n Aj) = 0. 
Thus, there is no SO. However, if we take two particles, 
we can choose x = o~ z ® o~ y and y — a y ® o~ z , so that the 
SOP is one. In general, we have Rl(u) = 1 as expected. 



2D Systems: Membrane order. Since the exis- 
tence of a local symmetry, and the possibility of deter- 
mining it numerically via the quantity Rl defined in (| 14[) . 
seems to be an appropriate definition of SO in ID sys- 
tems, we will now try to extend it to 2D. In this case 
there is more freedom in the choice of locations for the 
local unitaries u. First we can let them act on the whole 
lattice. If this leaves the state invariant (or the respective 
Rl is not vanishing) we will say that we have membrane 
order (MO). We can also put them as a string of opera- 
tors, in which case we will talk about SO [20J, or even in 
more sophisticated configurations, as a band of operators 
(BO). Note that the MO so defined shares all the desired 
properties for the SO in ID and thus provides a natu- 
ral generalization to higher dimensions. In particular, it 
should not exhibit spontaneous breakdown when switch- 
ing on couplings in the second dimension as pointed out 
in for the SO. 

To gain more inside we will consider the generaliza- 
tion of MPS /FCS called projected entangled pair states 
(PEPS) [H, [l3| where the matrices Ak are replaced by 
tensors Bk whose degree depends on the geometry of the 
lattice. For these states local symmetries can arise in a 
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a) 
b) 



V <= A^wV A^wV A^wV A=t>V 

<^ <?> 



V <= B=^>V 
V 



W 1 

+ V + 

V<=Bi=>V V^B^V 



w 



1 



FIG. 1: (a) Symmetries in FCS/MPS. The tensor A has three 
indices, one corresponding to the physical spin (pointing in) 
and two for the virtual systems (pointing left and right). Ap- 
plying U to the physical index is equivalent to applying V 
to the virtual ones. Since the tensor are contracted in a line 
the unitaries V and V' cancel and thus the state does not 
change, (b) 2D generalization (square lattice). The tensor 
B has one physical and four virtual indices. When applying 
U to the former, we can have different effects on the virtual 
system which will, after contraction, leave the state invariant. 
The example on the right allows for SO. 



similar way as in the ID case (FigfJ}. We interpret each 
tensor as an operator B — \s){(p s \ : (-ffo)® 4 ^ H$, 
where H$ (Hp) is the Hilbert space corresponding to 
the physical (virtual) spins. The PEPS exhibits a local 
symmetry if there exist unitaries V& and U ^ 1 with 



UB = BVi 



(16) 



such that when contracting the indices of B to create the 
state 'J (analogous to the matrices A in Eq.Q) the V's 
cancel [see Fig.Q]. Thus, we will have (y\U® N \^) = |*) 
and therefore MO. Note that we have the possibility of 
having different V's, in contrast to what happens in ID. 
This structure also allows us to understand the (dis-) 
appearance of SO or BO. However, in the 2D case, the 
connection between the existence of a local symmetry 
and Eq. (fT6|) is _less straight and will be analyzed in 
detail elsewhere 



14] 



Example 3: AKLT state [ll|. In this case 5 = 2, 
\<j> a ) = (<Jy ® <j y ® 1 <g> l)|0~ s ) an d 0s is an orthonormal 
basis of the symmetric subspace of (C 2 )® 4 . U and V 
correspond to a 5 and 2-dimensional representation of 
SU(2), and thus fulfill the condition of MO. One can 
prove that there is no SO even for a simple ladder formed 
by two chains. 

Example 4: Cluster state [l5|. Here, 5 = 1/2 and 

|0-i/ 2 ) = l + +oo), \4> l/2 ) = I — ii) with |±) = |o)±|i). 

One can take U = a x , V u = V r = <r x and Vd = Vi = a z . 
But one can also take U — a z , three of the V's equal 
to 1 and decide either V u = a z , V r — cr z , V = a x or 
Vi = a x . As a consequence, there is MO. One can readily 
show that there is no SO. However still there is a BO, in 
the sense that this state is an eigenstate of the operator 
obtained by applying unitaries to three consecutive lines 
(a z to the first and the third and a x to the middel one). 



Example 5: Toric code [16|. This is a state with 
5=1/2 and we have alternative tensors in the A and 
B sublattices, In this case we have alternating projec- 
tors B = |0)<$ + $ + |« Wr + |l)($_$_| uWr and B' = 
|0)($ + $+|„ rM + \l){®-3>-\ ur id, resp. One can take 
U = a x , Vi = V r = a z and V u = Vd = 1. Since there is 
no unitary to be cancelled in the up and down positions, 
it has SO (and also MO). 

Conclusion. We have shown that the existence of 
string order is, in the framework of FCS, equivalent to 
the existence of a local symmetry. This gives a direct 
explanation of many of its intriguing features, like its 
robustness/fragility under perturbations or the capabil- 
ity of detecting phase transitions. Though the FCS case 
provides evidence for the generality of this equivalence 
it remains an open problem to extend this beyond FCS, 
e.g., to all ground states of gapped local Hamiltonains. 
We have characterized the existence of local symmetries 
in FCS, related it to intertwining isometries and covari- 
ant channels (CI and C2), and sketched the case of 2D 
systems through PEPS. The list of examples can easily 
be extended and shows, for instance, in a simple way the 
presence of string order in topological ordered states, as 
it is illustrated in Kitaev's toric code. 

The present work is an example that MPS/FCS and 
PEPS are not only useful for numerical algorithms, but 
also to prove and clarify interesting statements. 

The obtained results shed light on the role of sym- 
metries in spin systems also in other contexts like 
Lieb-Schultz-Mattis type theorems [l7[. The respec- 
tive relation between integer/half-integer spins and ir- 
reducible/reducible representations will be discussed in 
a forthcoming paper 14| together with a more detailed 
investigation of the conditions for symmetry. 
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[1] M. den Nijs and K. Rommelse, Phys. Rev. B 40, 4709 
(1989). 

[2] E. G. Dalla Torre, E. Berg, E. Altman, Phys. Rev. Lett. 
97, 260401 (2006); F. Anfuso, A. Rosch Phys. Rev. B 75, 
144420 (2007). 

[3] Eugene H. Kim, G. Fath, J. Solyom, D. J. Scalapino 
Phys. Rev. B 62, 14965 (2000); S. Todo, M. Matsumoto, 
C. Yasuda, H. Takayama, Phys. Rev. B 64, 224412 
(2001). 

[4] F. Anfuso, A. Rosch Phys. Rev. B 76, 085124 (2007). 



5 



[5] D. Perez-Garcfa, F. Verstraete, M.M. Wolf, J.I. Cirac, 

Quant. Inf. Comp. 7, 401 (2007). 
[6] M. Fannes, B. Nachtergaele and R. F. Werner, Commun. 

Math. Phys. 144, 443-490 (1992). 
[7] G. Vidal, Phys. Rev. Lett. 91, 147902 (2003). 
[8] F. Verstraete, M.A. Martin-Delgado, J.I. Cirac, Phys. 

Rev. Lett. 92, 087201 (2004). 
[9] F. Verstraete, J.I. Cirac, Phys. Rev. B 73, 094423 (2006); 

M. B. Hastings J. Stat. Mech. (2007) P08024. 
[10] F. Verstraete, D. Porras, J.I. Cirac, Phys. Rev. Lett 93, 

227205 (2004). 

[11] I. Affleck, T. Kennedy, E. H. Lieb and H. Tasaki, Com- 
mun. Math. Phys. 115, 477 (1988). 

[12] F. Verstraete, J.I. Cirac, |arXiv:cond-mat /0407066| 

[13] D. Perez-Garcfa, F. Verstraete, M.M. Wolf, J.I. Cirac, 
larXiv:0707.2260l 

[14] In preparation. 

[15] R. Raussendorf, H.J. Briegel, Phys. Rev. Lett. 86, 5188 
(2001); F. Verstraete, J.I. Cirac, Phys. Rev. A 70, 



060302 (R) (2004). 
[16] A. Yu. Kitaev, Annals Phys. 303, 2 (2003). 
[17] E. Lieb, T. Schultz, and D. Mattis, Ann. Phys. 16 (1961), 
407466; M. B. Hastings, Phys.Rev. B 69 (2004) 104431; 
B. Nachtergaele, R. Sims, Commun. Math. Phys. 276 
(2007) 437-472. 
[18] In the context of MPS one finds an alternative but equiv- 
alent definition for string order [8J involving a finite chain 
with access to the virtual levels in the endpoints. 
[19] If one perturbs the Hamiltonian without keeping the sym- 
metry, it is clear from our picture that the string order 
will vanish, even if the gap does not close. This explains 
the fragility of the SO (cf. @]) 
[20] In this case one can obtain again the equivalence with the 
'classical' definition of string order. In order to see this, 
it is enough to invoke the injectivity condition defined in 




